function Hv = H( a, b )
% H Возвращает вектор фокусировки дл нашей антенной системы
%  Входные аргумента - азимут и угол места

[kx, ky, kz] = sph2cart(b, a, 1);
k= [kx; ky; kz];

r{1} = [1/4, -1/4, 0];
r{2} = [-1/4, -1/4, 0];
r{3} = [-1/4, 1/4, 0];
r{4} = [1/4, 1/4, 0];

NAnt = length(r);

phi = nan(NAnt, 1);

for i = 1:NAnt
 phi(i) = 2*pi*r{i}*k;
end

Hv = exp(1i*phi);
end

